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ABSTRACT 

We present fully general relativistic equations governing viscous transonic flows 
in vertical equilibrium in Kerr geometry. We find the complete set of global solutions 
(both for Optically thick and optically thin flows) in the weak viscosity limit. We show 
that for a large region of parameter space, centrifugal pressure supported standing 
shocks can form in accretion and winds very close to the black hole horizon, both for 
co-rotating and contra-rotating flows. We compute the nature of the shear tensor for 
complete transonic solutions and discuss the consequences of its reversal properties. 
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1 INTRODUCTION 

The physics of accretion processes on black holes has fas- 
cinated astrophysicists for last two decades (Novikov & 
Thorne 1973, hereafter NT73; Lynden-Bell, 1978; Moncrief, 
1980) but the full complexity of the flow behaviour is yet 
to be unfolded in a comprehensive manner. For example, 
one of the aspects that is recently studied more actively, 
is the accretion flow with 'discontinuities' or shock waves. 
The important numerical works of Wilson (1978) and Haw- 
ley, Smarr & Wilson (1984, 1985) show that sub-Keplerian 
flows experience centrifugal barrier close to the black hole 
and produce shock waves. Though these fully general rela- 
tivistic simulations revealed the existence of only traveling 
shock waves, subsequent theoretical and numerical study us- 
ing pseudo-Newtonian potentials show that shocks are in- 
deed stationary in a large region of the parameter space 
(Chakrabarti, 1989 [C89]; 1990a,b [C90a],[C90b]; Molteni, 
Lanzafame & Chakrabarti, 1994; Sponholz & Molteni, 
1994). Similarly, the formation of sub-Keplerian flows from 
Keplerian disks in radiation (Paczyhski & Wiita, 1980) and 
ion pressure (Rees et al., 1982) supported disks have been 
shown only in the context of pseudo-Newtonian geometry 
(C90ab, Chakrabarti, 1996 [C96]). So far, no global solu- 
tion of rotating, viscous, transonic flow (with cooling effects 
included) has been obtained using full general relativistic 
considerations. Only inviscid solutions of conical flows (with 
simplistic accretion rate equation) are present (C90b). 

In the present paper, we derive the equations governing 
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the viscous transonic flows in Kerr geometry and find global 
solutions using the usual sonic point analysis. We present ex- 
plicit examples of solutions in weak viscosity limit. We see 
that similar to our earlier studies using peudo-potentials, 
standing shocks can form around Kerr black holes for a 
large range of initial flow parameters. The maximum num- 
ber of sonic points remains three as was first pointed out 
by Liang and Thompson (1980) in the context of equatorial 
accretion in Schwarzschild geometry. The important differ- 
ence between the present study and the earlier works (e.g., 
Chakrabarti 1989) is that the present work does not use 
pseudo-potential and therefore clearly satisfies the bound- 
ary condition on the horizon. Secondly, shock locations, 
particularly in prograde flows, are very close to the black 
hole, roughly at half as much distance as those obtained in 
Schwarzschild geometry. In retrograde flows, the shocks form 
farther away. These may have important implications in the 
study of hard radiations from black hole candidates (e.g., 
Chakrabarti & Titarchuk, 1995). Also, for the first time, 
the nature of shear and angular momentum transport in a 
viscous transonic flow just outside the black hole horizon are 
discussed. 

In the next Section, we present the basic equations. 
In §3, we solve these equations in weak viscosity limit and 
show examples of global solutions with or without stationary 
shocks. In §4, we present the behavior of shear tensor which 
is responsible for angular momentum transport. Finally, in 
§5, we make concluding remarks. 
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2 BASIC EQUATIONS 



where, 



We use t, r, 9 and z as the coordinates. We choose the geo- 
metric units where G — Mbh = c = 1 (G is the gravitational 
constant, M bh is the mass of a black hole and c is the veloc- 
ity of light). We also consider \6 — tt/2\ << 1 for a thin flow 
in vertical equilibrium. We consider a perfect fluid with the 
stress-energy tensor, 



(1) 



where, p is the pressure and p = po(l+n) is the mass density, 
7r being the internal energy. We ignore the self-gravity of the 
flow. We assume the vacuum metric around a Kerr black hole 
to be of the form (e.g., NT73), 



ds — g llv dx IJ, dx v 



Where, 



■ 7 —^-dt 2 + -4r (dd> — iodt) 2 + T -rdr 2 + dz 2 
A r 2 A 

(2) 
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Here, g^ v is the metric coefficient and w M is the four velocity 
components: 

1/2 

(3a) 



u t 
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(l-V 2 )(l-M)(g^ + lg t4> ) 



and 



U4, = —lu t 
where, the angular velocity is 



g±4> + igu 

94><t> + 



(3b) 



(4) 



and I is the specific angular momentum. The radial velocity 
V in the rotating frame is (C90b, p. 137) 



V 



(I -ill) 1 / 2 



where, 



v = (- 



UrU 1/2 
UtU* 



Since V — 1 on the horizon and even in the extreme equation 
of state of p = p/3, the flow must be supersonic on the hori- 
zon. Thus, any black hole accretion is necessarily transonic 
(also, see C90b). 



2.1 Equations Governing a Viscous Transonic 
Flow 

In the present paper, we shall concentrate on the time in- 
dependent solutions of the underlying hydrodynamic equa- 
tions. The equation for the balance of the radial momentum 
is obtained from (w M u„ + g ilv )T^' — 0: 



dr r A r a r° r 1 P + P dr 





(5) 



7 2 = [l-^4(n-a,)]- 1 , 



B = (fia-1) 2 -fiV, 



and 



Here and hereafter we use a comma to denote an ordinary 
derivative and a semi-colon to denote a covariant deriva- 
tive. The baryon number conservation equation (continuity 
equation) is obtained from (po« M ) ;M = which is, 



where, 



M = 27rrtfE = 2nr-dp H 



H _ ( P V /2r 3/2 Ar 2 + a 2 f-Aa 2 1/2 
° { p ' 7 [ (r 2 + a 2 ) 2 + 2Ao 2j 



(6) 



is the height of the disk in vertical equilibrium (NT73). The 
equation of the conservation of angular momentum is ob- 
tained from {5%T tlv );v — 0, and one obtains, 



where, 



r) = vpo 



is the coefficient of dynamical viscosity and v is the coef- 
ficient of kinematic viscosity. When rotation is dominant 
(# < u ) the relevant shear tensor component crj is given 
by (Anderson & Lemos, 1988), 



A 3/2 7 3 f2. r A 1/2 
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(7) 



so that the angular momentum equation takes the form, 

£-£ + = -l^W 2 . (8) 

Here we have modified earlier work of NT73 and Peitz (1994) 
in that the fluid angular momentum C = —hu^, rather 
than particle angular momentum — is used. That way, 
for an inviscid flow (77 = 0) one recovers C — constant 
as in a fluid picture. Similarly, the radial velocity term 
is included (eq. 3) and angular momentum is allowed to 
be non-Keplerian (eq. 8). £+ is the angular momentum 
on the horizon since the rotational shear (as defined by 
cq. 7) vanishes there. In presence of significant radial ve- 
locity, the shear in eq. (7) is to be replaced by its full 
expression, a"" = (u^P (iv + u^P^)/2 - QP^ /3 where 
P M " = g^ v + u^vl" is the projection tensor and O = is 
the expansion (e.g., Shapiro & Teukolsky, 1983) which no 
longer vanishes on the horizon (see, Fig. 6 below). In that 
case, C + will no longer be the specific angular momentum of 
the flow at the horizon, instead it will be at the place where 
the net shear effect vanishes. In Section 4, we shall compare 
these shear components in a realistic solution. 

Entropy generation equation is obtained from the first 
law of thermodynamics along with the baryon conservation 
equation (5") ;M = [2-qa^a^yT - Q~ : 



dh , 
dr po dr 



= 2uT,au V cr' J 
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where Q + and Q~ are the heat generation rate and the 
heat loss rate respectively, h is the specific enthalpy: h = 
(p + p)/po- Here, we ignore the terms contributed by radia- 
tion. Using rotational shear as given in eq. (7) , the entropy 
equation takes the form, 



dh 
dr 



po dr 



-Q- 



(9) 



Of course, for accuracy, we use the full expression for crj as 
discussed above. 

This set of equations are solved simultaneously keeping 
in mind that the shock waves may form in the flow, where, 
the following momentum balance condition 

W-n+(W-+Y, -){u> 1 _n IJb )u_ = IF + n l/ + (W/ + +E 0+ )(<n, 1 >+ 

(10) 

along with the conservation of energy and mass fluxes (to- 
gether, these conditions are known as the Rankine-Hugoniot 
conditions) must be fulfilled. Here, n M is the four normal vec- 
tor component across the shock, and W and E are vertically 
integrated pressure and density on the shock surface. The 
subscripts — and + denote the pre-shock and post-shock 
quantities respectively. 

The equations presented above are applicable to opti- 
cally thin as well as optically thick flows for any general heat- 
ing and cooling processes. For a given viscosity prescription 
and the exact cooling processes (depending on the optical 
depth of the flow), it is usual to reduce the above set of 
equations in the form (Bondi, 1952): 



du_N_ 
dr ~ D 



(11) 



where TV and D are the smooth functions of radial coordinate 
(unless there are non-linearities which prevent such reduc- 
tion. In that case sonic curve analysis is done, see Flammang, 
1982). The procedure of obtaining the complete solution is 
then similar to what is presented in obtaining the global 
solutions of viscous transonic flow using pseudo-Newtonian 
potential (C90a,b; C96 and references therein). Presently, 
we concentrate on the solution of the above set of equations 
in the inviscid limit. Full discussion on strongly viscous flows 
will be presented elsewhere. 



2.2 Sonic Point Conditions in Weak Viscosity 
Limit 

In the case when viscosity is negligible, equations (5) and 
(8) could be integrated to obtain the energy and angular 
momentum conservation equations (Moncrief, 1980; Flam- 
mang, 1982; C90b), 



and 



7 P + P 

hu t — u t 

Pa 



p±P i 
Po 



-c. 



(12a) 



(126) 



From eq. (6) we get the constant accretion rate as, 

M = 2-Krtfp H . (12c) 

In the weak viscosity limit, the entropy equation (9) can 
be replaced by an adiabatic equation of state P = Kp\ +1 ^ n . 
K, which is a measure of entropy, is constant in the flow 



but it changes at the shocks (the entropy generated may be 
either radiated or advected away depending on the cooling 
effciency of the flow at the shock; see C90b) , n is the poly- 
tropic index related to the adiabatic index F by F = 1 + 1/n, 
and is assumed to be constant throughout the flow. The spe- 
cific enthalpy and the density in terms of the sound speed 
a 3 = (dP/dp) 1 ^ 2 become, 



h - 



Po = 



K(n + l)(l-na 2 )' 



(13a) 



(136) 



First we rewrite the energy and mass conservation equa- 
tions on the equatorial plane as (9 — ir/2), 



1 /AX 1 / 2 



= ^ I ( 

1-na 2 {l-V 2 ) 1 ' 2 \V) 



(14a) 



and 



as vn+1/2 V f Ar 3 {(r 2 + a 2 ) 2 - Aa 2 } 
1-na 2 ' (1 - V 2 ) 1 / 2 [ 7 {(r 2 + a 2 ) 2 + 2Aa 2 } 



1/2 



where, 



(146) 



2Z 2 a, 



V = (1 - m){g u + Ifto) = r 2 + — (1 - j) 2 - l 2 {l - — 



(15a) 
(156) 



The quantity M will be called the entropy accretion rate 
as Ai oc K n M. It usually assumes two different values at 
two saddle type sonic points and the one with the smallest 
M joins horizon with infinity (C89, C90b). This quantity, 
which is very useful in the study of shocks, was defined for 
the first time in C89 and was simply termed as 'accretion 
rate'. In presence of an accretion shock, the transonic flow 
passes through both the sonic points. 

Differentiating equations (14a) and (14b), and eliminat- 
ing terms involving derivatives of a s , we obtain, 



dV N cia 2 — C2ai 
dr D bia2 — &2ai 



where, 



ai = 



a 2 



62 



2na s 
1 — na 2 

2n+ 1 
a s (l - na 2 ) 

V 

~ 1 - V 2 
1 



ci = 



V(l - V 2 ) 
r-1 r-^(l-f) 2 



3 r 

C2 = — + - 

2r 



V 



A 7 #0 



(16) 

(17a) 
(176) 
(17c) 
(17d) 

(17c) 
(17/) 
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Where a prime (') denotes ordinary derivative with respect 
to r. From the vanishing condition of N and D at the sonic 
points, one obtains the so called sonic point conditions as, 

C2 

and 



(18a) 



= V C ( 



2n + 1 .1/2 
2n ' 



(186) 



Here, the subscript 'c' refers to the quantities evaluated at 
the sonic point r = c. The sonic point condition (7b) im- 
plies a Mach number M = V/a s = (2n/2n + 1) 1/2 at the 
sonic point. This is of the same form as the flow in vertical 
equilibrium in pseudo- Newtonian potential (C89), though 
the velocities here are defined with relativistic corrections. 
In conical flows in Kerr geometry (with a simpler accretion 
rate equation), M = 1 at the sonic point (C90b). 

In obtaining a global solution one supplies the conserved 
quantities at the inner (such as the disk surface) or the outer 
(e.g, Keplerian or sub-Keplerian flows injected at the outer 
region) boundary, depending upon whether one is interested 
in the wind solution or the accretion solution. For a given 
angular momentum I, the remaining unknowns are V(r), 
a 3 (r). But one requires only one extra boundary condition, 
e.g., £, since two sonic point conditions (eqs. 18a and 18b) 
introduce only one extra unknown, namely, r c . Thus, the 
supply of the initial specific energy £ and the specific angu- 
lar momentum I are sufficient for a complete solution from 
the horizon to infinity. For a viscous flow, one clearly has to 
supply the distribution of viscosity rj(r) (e.g., ion or mag- 
netic viscosity) itself. Simple Shakura & Sunyaev (1973) vis- 
cosity prescription may not be very useful since that stress 
—aP is always negative, while in a general relativistic flow 
stress can change signs. Note that by definition, £1 = u> on 
the horizon and thus the flow co-rotates with the black hole. 
Instead of specifying various quantities at the flow bound- 
ary, one can alternatively specify the location of a critical 
point along with the energy or the angular momentum as in 
C90a,b. 

In much of the parameter space, the flow is expected 
to be smooth as in a Bondi flow. However, if the angular 
momentum is significant, matter can pile up behind the cen- 
trifugal barrier close to the black hole and form a standing 
shock wave. At the shock, apart from the continuity of the 
energy and mass flux, the relativistic momentum balance 
condition (eq. 10) must be satisfied. Using Newtonian defi- 
nition of the vertical integration (since we are dealing with 
thin flows anyway) as in C89 and the definition of entropy 
accretion rate M, we find that at the shock, the following 
quantity: 



n = 



H3/2 



( 2 I V 

\3T-1 + ^JTT 



V 2 ) 



M 



(19) 



should be continuous. One can use different prescriptions, 
but the general behaviour is not expected to change. 



3 BEHAVIOUR OF SONIC POINTS AND 
GLOBAL SOLUTIONS 

Before obtaining a globally complete solution, we can make 
a general study of the sonic point behaviour of a weakly vis- 



cous flow. Using eq. 14a, eq. 18a and eq. 18b, we can compute 
the specific energy £ of the flow (which is conserved) as a 
function of the sonic radius r c and specific angular momen- 
tum of the flow. Figs, l(a-b) shows our findings for a = 0.95. 
We choose T = 4/3, i.e., n = 3. In Fig. la, from the top curve 
to the bottom we use, I = 2.0, 2.1, 2.2, 2.3, 2.4 and 2.5. 
The number of sonic points of the flow for a given energy 
and angular momentum will depend on the number of inter- 
sections between a £ = constant line and a constant angular 
momentum curve. If the slope of the curve is negative the 
sonic point is saddle type and is physical, otherwise it is of 
circle type or 'O' type (see, C89, C90a,b) which is unphysical 
and irrelevant. We note that for a very high energy £ » 1 
or for low angular momentum there is only one saddle type 
sonic point as in Bondi accretion (Bondi, 1952). But near 
£ — 1 (i.e., the rest mass), there could be three sonic points 
(see Fig. 3 below). The innermost and the outermost inter- 
sections represent saddle type sonic points and the middle 
intersection represents the 'O' type sonic point. For low an- 
gular momentum (say, for / = 2.0) there is only one sonic 
point for all energy with which matter enters the black hole. 
This is similar to a spherical Bondi flow. However, for higher 
angular momentum, say I = 2.3, there could be two saddle 
type sonic points for a range of energy. The presence of the 
second saddle type point is the general relativistic effect, as 
discussed in C90b or C96. In Fig. lb, we present results for 
retrograde flows with angular momenta (from the top curve 
to the bottom), / = 3.4, 3.8, 4.2, 4.6, 5.0 and 5.4. Here, 
to have three sonic points the flow requires a larger angular 
momentum. This is because the matter has to fight against 
the frame-dragging effects which forces the flow to co-rotate 
with the horizon. The terms containing a/l in Eq. (15a) and 
Eq. (17e) represent the coupling of the spin angular momen- 
tum of the black hole and the orbital angular momentum of 
the particle which bring about this change of behaviour. In 
either cases, however, the flow angular momentum that we 
consider here are around marginally stable value or even 
less. Thus, they are not high in any real sense. 

For a higher polytropic index, say, for T = 5/3, the 
number of saddle type sonic points is only one. The T — 
Tcrit at which the transition from two saddle type points to 
just one saddle type sonic point takes place depends on the 
geometry of the flow. Typically, we find that for a flow in 
vertical equilibrium, Y crit ~ 1.5 (see also, C90b; C96). 

Another important property of the flow which could 
be studied before solving the full set of equations is the 
dependence of entropy accretion rate on specific energy at 
various sonic points for a given angular momentum of the 
flow. Since one requires only two free parameters for a global 
solution, when a pair (£, I) is supplied at the outer boundary, 
M must be regarded as an eigenvalue of the problem. Note 
that only the product of entropy function K n and the mass 
accretion rate M is determined by the boundary condition, 
not the individual quantities. Thus, a specific flow property, 
such as shock location, could be completely independent of 
the mass flux M. When cooling processes (which depend 
on mass accretion rate) are used, the result will depend on 
accretion rates, except when cooling relation is simple (such 
as Q~ /Q + ^constant, see C96). 

In Fig. 2, we present the variation of the specific en- 
ergy £ with entropy accretion rate M as the sonic point is 
changed. This is for inviscid flow (r)(r) = 0). The direction 
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Fig. la: Plot of specific energy £ as functions of the sonic point 
(x c ) for I = 2.0, 2.1, 2.2, 2.3, 2.4 and 2.5 (from top to bottom). 
a = 0.95 is used. For some range of energy, the curves of certain 
I may have three intersections, signifying the presence of three 
sonic points and possibility of shock waves. 
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0.85 - 



0.8 1 1 1 1 1 1 1 1 1 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 
Log (Sonic Radius) 

Fig. lb: Same as (a) but for contra-rotating flow with I = 
3.4, 3.8, 4.2, 4.6, 5.0 and 5.4 (from top to bottom). The sonic 
points are located farther out in this case. 

of the arrow on the curve indicates the direction in which 
the sonic point is increased. The specific angular momen- 
tum of the flow is I = 2.3 and a — 0.95 is chosen as before. 
The branches AMB and CMD are drawn for the inner and 
the outer sonic points respectively and the branch BC is 
drawn for the circle type sonic point (C89). For a given mass 
accretion rate M, the entropy accretion rate measures the 
specific entropy of the flow. In the energy range spanned by 
MC, Mo > Mi, and in the energy range spanned by MD, 
Mo < Mi- (Here, the subscripts i and o denote the quan- 
tities at the inner and the outer sonic points respectively.) 
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0.85 '— 1 1 1 1 

0.0001 0.0002 0.0003 0.0004 

Entropy Accretion Rate 

Fig. 2: Plot of specific energy £ against the entropy accretion 
rate M at all the sonic points of the flow of angular momentum 
I = 2.3 around a Kerr black hole of a = 0.95. Sections AMB 
and CMD denote results on inner and outer saddle type sonic 
points respectively, and the section BC denotes results of center 
type sonic point. Horizonal arrows aa, and ww schematically show 
shock transitions in accretion and winds respectively. Arrows on 
the curves show the direction in which the sonic point is increased. 

However, since only smaller of the two Ms joins the flow 
at infinity (C90b), only allowed (shock free) flow must have 
parameters either on AM (accretion or wind through inner 
sonic point) or on MD (accretion or wind through outer 
sonic point). In some regions of the parameter space close 
to M, two different transonic flows (e.g., two flows with the 
same energy and angular momentum, but different entropy 
accretion rate) can be 'glued' together through shock con- 
ditions. For example, a flow, first passing through an outer 
sonic point on MD can undergo a shock aa and generate 
enough entropy to 'land on' the branch MB so that it can 
now pass through the inner sonic point. This is an exam- 
ple of the accretion shock. We shall show this solution in 
Fig. 3 below. Conversely, a wind flow can first pass through 
an inner sonic point on MA and after passing through a 
shock can land on a point on the branch MC so that it can 
pass through the outer sonic point also. This solution is also 
shown in Fig. 3 below. In order to have a shock, of course, n 
of Eq. (9) must be continuous across the shock (that defines 
'some' region of energy close to the point M on the curve.). 

So far, we discussed the nature of solutions when both 
the saddle type sonic points are present. This is true only 
in an energy range AS spanned by the curve CMD (see, 
Fig. 3 below). Otherwise, the flow will have only one sad- 
dle type sonic point, and it will surely not have a standing 
shock wave of the kind we are discussing here (though shocks 
by external irradiation which forces the flow to be subsonic 
cannot be ruled out). The whole range of AS does not allow 
shock solution, but only a part, in which Rankine-Hugoniot 
condition is satisfied, is allowed. This will be shown below. 

In Fig. 3 we classified the entire parameter space ac- 
cording to the type of inviscid solutions that is prevalent. 
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Fig. 3: Classification of the parameter space (central box) in 
the energy-angular momentum plane in terms of various topol- 
ogy of the black hole accretion. Eight surrounding boxes show 
the solutions from each of the independent regions of the param- 
eter space. Each small box shows Mach number M against the 
logarithmic radial distance r. Contours arc of constant entropy 
accretion rate M. Parameters from region N does not form any 
transonic solution. See text for details. 

We choose a — 0.5. (For classification of flows in pseudo- 
Newtonian geometry, see C89, C90ab). The adiabatic index 
r = 4/3 has been chosen. In the central box, we divide the 
parameter space spanned by (I, £) into nine regions marked 
by N, O, NSA, SA, SW, NSW, I, O* , I*. The horizon- 
tal line at £ = 1 corresponds to the rest mass of the flow. 
Surrounding this parameter space, we plot various solutions 
(Mach number M — v r /a s vs. logarithmic radial distance 
where v r is the radial velocity and a s is the sound speed) 
marked with the same notations (except N). Each of these 
solution topologies has been drawn using flow parameters 
from the respective region of the central box. The accretion 
solutions have inward pointing arrows and the wind solu- 
tions have outward pointing arrows. The crossing points are 
'X' type or saddle type sonic points and the contours of cir- 
cular topology are around '0' type sonic points. If there are 
two 'X' type sonic points, the inner one is called the inner 
sonic point and the outer one is called the outer sonic point. 
If there is only one 'X' type sonic point in the entire solution, 
then the terminology of inner or outer is used according to 
whether the sonic point is closer to or farther away from the 
black hole. The solutions from the region 'O' has only the 
outer sonic point. The solutions from the regions NSA and 
SA have two 'X' type sonic points with the entropy density 

50 at the outer sonic point less than the entropy density 

51 at the inner sonic point. However, flows from SA pass 
through a standing shock since the Rankine-Hugoniot con- 
dition is satisfied. The entropy generated at the shock Si — S 
is advected towards the black hole to enable the flow to pass 
through the inner sonic point. Rankine-Hugoniot condition 



is not satisfied for flows from the region NSA. Numerical 
simulations show (Ryu, Chakrabarti & Molteni, 1997) that 
the flow from this region is very unstable and exhibit peri- 
odic changes in emission properties as the flow constantly 
tries to form the shock wave, but fails to do so. The so- 
lutions from the region SW and NSW are very similar to 
those from SA and NSA. However, S > Si in these cases. 
Shocks can form only in winds from the region SW. The 
shock condition is not satisfied in winds from the region 
NSW. This may make the NSW flow unstable as well. A 
flow from region I only has the inner sonic point and thus 
can form shocks (which require the presence of two saddle 
type sonic points) only if the inflow is already supersonic due 
to some other physical processes. Each solution from regions 
I* and O* has two sonic points (one 'X' and one 'O') only 
and neither produces complete and global solution. The re- 
gion 7* has an inner sonic point but the solution does not 
extend subsonically to a large distance. The region O* has 
an outer sonic point, but the solution does not extend su- 
personically to the horizon! When a significant viscosity is 
added, the closed topology of I* opens up as described in 
C90ab and C96, and the flow may join with a cool Keplerian 
disk with £ < 1. These special solutions of viscous transonic 
flows should not have centrifugally supported shock waves as 
they have only one inner sonic point. However, hot flows de- 
viating from a Keplerian disk or, sub-Keplerian winds from 
companions or, cool flows subsequently energized by mag- 
netic flares (for instance) will have £ > 1 and thus could 
have standing or periodically varying shock waves as dis- 
cussed above. The post-shock flow radiates most of the ob- 
served hard radiation as shown by Chakrabarti & Titarchuk 
(1995). Note that parameters from region N does not pro- 
duce any transonic solution. 

In C89 and C90ab, it was found that shock conditions 
were satisfied at four locations: r s i, r S 2, r S 3, r S 4, though r s i 
and r 3 4 were found to be not useful for accretion on black 
holes. Out of r S 2 and r S 3, it was shown that r s s is stable 
for accretion flow and r S 2 is stable for winds (Chakrabarti 
& Molteni, 1993, also see, Nakayama, 1994 and Nobuta & 
Hanawa, 1994). We plot here only r s s here in SA and r S 2 in 
SW solutions. Here, o and i are the outer and inner sonic 
points respectively. In the box containing a solution from 
SA, we chose a = 0.5, I = 3, £ = 1.003. For these param- 
eters, the eigenvalue of the critical entropy accretion rates 
at the two saddle type sonic points are Mi = 2.74 x 10~ 5 
and Mo = 1.491 x 10~ 5 respectively. Here, M a < Mi, 
hence the flow through the outer sonic point joins the hori- 
zon with infinity (single arrowed curve). The flow forms a 
shock and jump onto the branch which passes through i 
as shown by double arrows. The stable shock (shown by 
a vertical dashed line) is located at 03 = r S 3 — 32.29 (in 
C89 notation). Only this jump, namely, a generation of en- 
tropy of amount Mi — Mo is allowed in order that the tran- 
sonicity of the post-shock flow is guaranteed. The entropy 
generated at the shock is advected through the inner sonic 
point. The flow is inefficiently cooled, which keeps the en- 
ergy of the flow constant. This makes the flow much hotter 
than a Keplerian disk. (This is typical of advective disks. 
See, C89, C96.) In SW solution, we chose a = 0.5, / = 3, 
£ = 1.007. For these parameters, the eigenvalue of the crit- 
ical entropy accretion rates at the two saddle type sonic 
points are Mi = 3.12 x 10~ 5 and M = 5.001 x 10~ 5 re- 
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Fig. 4: Variation of shock locations with energy in accretion and 
winds for various specific angular momenta I (marked on curves). 
a = 0.5 is chosen. Segments marked 03 and W2 (solid curves) 
represent stable shocks in accretion and winds respectively. Other 
two segments (02 and W3) represent formal shock locations which 
are unstable. 

spectively. Here, Mi < Mo, hence the flow through the 
inner sonic point i joins the horizon with infinity (single ar- 
rowed curve). The accretion flow can no longer form a shock. 
But a wind, first passing through i can, as shown in double 
arrows. The stable shock (shown by a vertical dashed line) 
is located at W2 = r a 2 = 6.89 in this case. Only this jump, 
namely, a generation of entropy of amount Mo — Mi is al- 
lowed at the shock in order that it can escape to infinity 
through the outer sonic point O. This consideration, along 
with the continuity of II (Eq. 19) allows one to locate sta- 
tionary shock waves in a flow. Note that though the flow has 
a shock-free solution (passing through o for accretion in SA 
solution and through i for winds in SW solution in Fig. 3), 
the flow would choose to pass through a shock because the 
latter solution is of higher entropy. This fact has been ver- 
ified through numerical simulations of accretion and wind 
flows (Chakrabarti & Molteni, 1993; Molteni, Lanzafame & 
Chakrabarti, 1994). It is to be noted that the angular mo- 
menta associated with solutions which include shocks arc 
not arbitrarily large. Rather, they are typically less than 
the marginally stable value Z ms as indicated in Fig. 3. 

Global solutions which contain shock waves are not iso- 
lated solutions, but are present in a large range of energy 
and angular momentum. In Fig. 4, we show the variation 
of shock locations as a function of specific energy £. Each 
set of curves, drawn for various specific angular momen- 
tum (marked on the set), consists of four segments: two 
for accretion (112 = r a 2 and 03 = r S 3) and two for winds 
(w2 = r S 2 and W3 = r S 3). As discussed above, 0,2 and W3 
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Fig. 5b: Example of shock transition in a retrograde flow. The 
parameters are: a = —0.95, I = 4.0, £ = 1.003. 

(dotted curves) are unstable while 0,3 and W2 (solid curves) 
are stable. Kerr parameter a — 0.5 is chosen. This exam- 
ple shows that stable shocks can form for a very wide class 
of flows. For corotating flows, the marginally stable and 
marginally bound angular momenta are l ms = 2.9029 and 
Imb = 3.4142 respectively. Thus the shocks form for angular 
momentum around these values. Since centrifugal barrier 
becomes stronger with angular momentum, shocks are lo- 
cated at larger radii for higher angular momenta. Another 
important point to note is that the shock location increases 
when the specific energy is increased. In a quasi-spherical 
flow, with the same input radial velocity and angular mo- 
mentum, the potential energy decreases with height (since 
the gravity becomes weaker), thereby increasing the spe- 
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cific energy. Thus, shocks in a three-dimensional flow arc 
expected to bend outwards. Such shocks have indeed been 
seen in numerical simulations of Hawley et al. (1984, 1985) 
and Molteni et al., 1996. Note that the parameters for stable 
shocks in accretion and winds are mutually exclusive. This is 
a consequence of the second law of thermodynamics, which 
requires that M + > M- (C89). 

We now compare the nature of shock solutions in pro- 
grade and retrograde flows. In Fig. 5a, we show a shock 
solution in prograde accretion with a = 0.95, I = 2.3, £ = 
1.001. The eigenvalues of the entropy accretion rates at the 
outer (O) and inner (/) sonic points are, Mo = 2.973 x 
10~ 06 , Mi = 4.108 x 10~ 5 respectively and the stable 
shock is located at r s z = 18.287. In Fig. 5b, we show a 
shock solution in retrograde accretion with a — —0.95, / = 
4.0, £ = 1.0025. The eigenvalues of the entropy accre- 
tion rates at the outer (O) and inner (I) sonic points are, 
Mo = 1.097 x 10~ 05 , Mi = 1.478 x 10" 05 and r s3 = 53.41. 
For a retrograde flow, shock locations are generally higher 
than that for a prograde flow. Not shown in the figure is an 
interesting fact that the Mach number of the flow remains 
perfectly finite (typically in the range 5 — 8 in our examples) 
on the horizon. 



4 GLOBAL NATURE OF SHEAR TENSOR 

In presence of viscosity, angular momentum transport is de- 
termined by the shear tensor. Whereas the computation of 
actual shear must be done in conjunction with the global 
solution of velocity components, presently we discuss its na- 
ture when the viscosity is weak, so that our inviscid global 
solution may be used instead. Around a Newtonian star, 
shear tensor component due to pure rotational motion can 
be defined as a r <t> = —rdQ/dr which is always positive 
(~ 3/2Q for a Keplerian disk), and hence the viscous stress 
is always negative (that is why —aP viscosity prescription of 
Shakura-Sunyacv [1973] was admissible). The viscous force 
in the <f> direction, caused by rubbing the adjacent layers 
generates a torque which would carry angular momentum 
only outwards. However, Anderson & Lemos (1988) pointed 
out that the shear can change sign near the horizon and 
they demonstrated this using velocity profiles of a cold ra- 
dial flow below marginally stable orbit. Below we show that 
this conclusion remains similar for our global solution. It is 
to be noted that this shear is 'passive', namely, computed 
after the velocity profiles are known from our inviscid so- 
lution. On the other hand, this is certainly more accurate 
than what is obtained by rotational velocity alone. In pres- 
ence of small viscosity, the behavior of our shear is expected 
to be the same. Computations with higher viscosities will be 
presented elsewhere. 

In Section 2, we presented two definitions of shear ten- 
sor, one for the rotationally dominated flow and the other for 
the general flow. In Fig. 6, we compare these two definitions 
(&$\rot and crj±) using the velocity profiles from our solu- 
tion. We also plot the gradient of angular velocity (dotted 
curve). The upper and lower panels are drawn using the 
parameters of Figs. 5a and 5b respectively. Thus they rep- 
resent typical behaviours in prograde and retrograde flows 
respectively. We plot only the solution branches which pass 
through the outer sonic point (solid curve represents for 



the supersonic branch, and the long dashed curve represents 
a r t j > _ for the subsonic branch). The other branches which pass 
through the inner sonic point behave in a similar manner. 
Also shock-free global solutions show a similar behaviour. 
The difference in these two cases is due to significant ra- 
dial velocity of the flow. We note that while the shear stress 
derived from rotational motion and that for the subsonic 
branch derived from general motion (crj_) vanish on the 
horizon, the general shear on supersonic branch (cr^ + , which 
is relevant for our work) does not vanish on the horizon. This 
shows that a flow just outside the horizon will transport 
angular momentum towards the horizon. Although for any 
realistic viscosity, this behaviour is not expected to make 
any difference in the angular momentum distribution, and 
therefore our earlier computations (e.g. C90ab) which ig- 
nore such effect (by assuming aP viscosity which is always 
positive, say) should remain valid. The reversal of the vis- 
cous stress is clearly due to the coupling of the rotational 
energy and the gravitational energy of the flow. Here, the 
'mass' equivalent of the rotational energy is also attracted 
by gravity, making the attractive force stronger than a that 
of a Newtonian star while the centrifugal force remains the 
same. The gravity always wins (giving rise to the well known 
'pit in the potential' for any angular momentum). To remain 
in equilibrium the flow angular momentum must also rise as 
it approaches the black hole horizon and to achieve this the 
angular momentum must be transported inwards. This is 
the origin of the trend of the shear tensor. The non-linearity 
(coupling among all the energy terms) which gives rise to 
this interesting behaviour is the hall-mark of general rela- 
tivity as is discussed in detail in Chakrabarti (1993). 



5 CONCLUDING REMARKS 

In this paper we presented the equations governing a vis- 
cous transonic flow in vertical equilibrium around a Kerr 
black hole and obtained the complete set of global solu- 
tions of them. We showed that for a large region of the 
parameter space, a transonic flow can have a shock wave as 
well. We showed that standing shocks can form much closer 
to the black hole than shocks in the Schwarzschild geom- 
etry. This causes the post-shock region to be much hotter 
and may thus contribute to very high energy X-rays and 
7-rays observed from black hole candidates (Chakrabarti & 
Titarchuk, 1995). Furthermore, the angular momentum re- 
quired to have shock waves are also lower, (typically I ~ 
GM bh /c for a ~ 1) than that required for a Schwarzschild 
black hole (typically / ~ 4GMbh/c for a = 0). This is sig- 
nificant, as even with a small angular momentum the flow 
can have a stationary shock. It is true that the region of 
parameter space which forms a shock is limited, but numer- 
ical simulations (Molteni, Lanzafame & Chakrabarti, 1994; 
Molteni, Ryu & Chakrabarti, 1996) indicate that flows even 
outside this range of parameter contain shock waves, partic- 
ularly because of the presence of turbulent pressure whose 
effects we have not included here. 

Earlier (C89, C90ab, C96) our studies were made using 
pseudo-Newtonian potential, where the the inner boundary, 
namely the horizon was not satisfactorily described. For in- 
stance, the angular velocity Q = l/r 2 was not equal to the 
angular velocity on the horizon even for a Schwarzschild ge- 
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Fig. 6: Comparison of rotational shear stress o r ^\rot (short dashed 
curves) with complete shear stress <r^ + (solid) along the super- 
sonic branch (passing through outer sonic point) for (a) a pro- 
grade flow (upper panel) and (b) a retrograde flow (lower panel). 
Also shown is d£l/dr (dotted curves). For comparison, results of 
the subsonic branch (<r£_) which passes through the outer sonic 
point is also shown (long dashed). Note the change in sign of shear 



near the horizon, 
and c r ± do. 



docs not vanish on the horizon, but a r . 



j ^_ uu. Here, subscripts '+' 
subsonic branches respectively. 



and '-' denote supersonic and 



ometry. Usually this was not a problem, since the potential 
energy was singular and the rotational energy due to any fi- 
nite angular momentum was negligible compared to it. How- 
ever, in our present analysis, such things do not arise. By 
definition, Q, = ju 1 co-rotates with the horizon and our 
solution continuously or discontinuously extends from the 
horizon to infinity. 

In this paper, we concentrated on the cases where the 
polytropic index n = 3, i.e., F = 4/3. When the flow is com- 
pletely gas pressure dominated, the flow may have r ~ 5/3 
in absence of cooling effects. In this extreme case, the flow 
will have only one (inner) saddle type sonic point (C90b) 
and the shocks can form only if the flow is already supersonic 
(e.g., when some wind is accreted). In reality, the flow will 
have intermediate polytropic index, dictated by cooling and 
heating processes and the possibility of shock formation will 
still be significant. The separation of topologies take place 
at around T — 1.5 for flows in vertical equilibrium (C96). 

We have discussed the behaviour of shear tensor compo- 
nents for our transonic solutions when the viscosity is negli- 
gible. Viscous stress shows a similar trend (namely, develop- 
ing a maximum) as in a cold, radial flow (Anderson & Lemos 
1988). However, unlike this case where the reversal occurs 
for a slowly rotating black hole, our hot, non-radial flow 
shows the trend of reversal even when the Kerr parameter 



is high. The conclusions we drew for an inviscid flow cannot 
change dramatically when a weak viscosity is present, ex- 
cept that instead of 'X' and 'O'. type sonic points, one will 
have 'X' and 'spiral' or 'nodal' type points (C90ab, C96). 
A flow which is cool and strictly Keplerian must be bound 
(8 < 1) and therefore will not have shocks (Fig. 3) with- 
out additional energizing mechanism. The Keplerian disk 
will simply become sub-Keplerian due to pressure effects 
and enter through the inner sonic point before disappear- 
ing into the black hole. On the other hand, if the Keplerian 
disk is hot, and/or the sub-Keplerian flow originated from a 
cool Keplerian disk is re-energized (by magnetic flaring and 
other effects, for instance), or, if the sub-Keplerian flows 
are mixed with supersonic winds coming from companions 
then they can have standing shocks provided viscosity is 
below a critical value (C96). In fact, even in an axisym- 
metric fully three dimensional flow, in full Kerr geometry, 
the shocks form more easily above and below the equatorial 
plane (Chakrabarti, in preparation) where the accretion is 
presumably highly advective and sub-Keplerian. These con- 
siderations were used by Chakrabarti & Titarchuk (1995) 
to prove that the transition of spectral states of black hole 
candidates is due to variation of accretion rates. The anal- 
ysis of 1100 days of BATSE data of Cyg X-l (Crary et 
al, 1996) shows an excellent agreement with such consid- 
erations. In the absence of shocks, the spectral properties 
may vary similarly (our conclusions regarding the produc- 
tion of weak hard tail component in converging flows are 
independent of whether shocks form or not) since viscous 
effects (or, equivalently, accretion rates) change the tran- 
sition radius where Keplerian disk becomes sub-Keplerian 
(C90ab; Chakrabarti & Titarchuk, 1995; C96). However, in 
this case variations in the hard and soft components are al- 
ways correlated which are not observed in most of the black 
hole candidates, especially in wind fed systems. Thus, we 
believe two independent accretion components, namely, Ke- 
plerian and sub-Keplerian, should be present rather than 
one single component becoming sub-Keplerian and produc- 
ing hard radiations. The most general solutions of advective 
disks in pseudo-Newtonian potential discussed these prop- 
erties (C96) and we showed in the present paper that the 
general conclusions do not change even when the computa- 
tions are carried out in full general relativity. 

In the present paper, we have written the most general 
equations in Kerr geometry which govern hydrodynamics of 
the flow (originally attempted by NT73 in the Keplerian 
limit) by including radial velocity, and correcting the angu- 
lar momentum transport equations with the use of conserved 
proper angular momentum and proper shear components. 
We already showed that the relevant shear components be- 
have very differently from that of a flow around a Newtonian 
star. In future, we will present solutions of these equations 
for highly viscous flows. 

The author thanks Kip Thorne and J. Peitz for discus- 
sion. 
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